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Heavy-meson semileptonic decays calculations on the lattice are reviewed. The focus is upon obtaining reliable 
matrix elements. Errors that depend upon the lattice spacing, a, are an important source of systematic error. Full 
0(o) improvement of matrix elements for arbitrary-mass four-component quarks is discussed. With improvement, 
bottom-quark matrix elements can be calculated directly using current lattices. Momentum dependent errors for 
0(a)-improved quarks and statistical noise limit momenta to around IGeV/c with current lattices. Hence, 
maximum recoil momenta can be reached for D decays while only a fraction of the maximum recoil momentum 
can be reliably studied for the light-meson decay modes of the B. Differential decay rates and partial widths are 
phenomenologically important quantities in B decays that can be reliably determined with present lattices. 



1. INTRODUCTION 



Semileptonic decays provide much needed in- 
formation about CKM matrix elements involv- 
ing heavy quarks. Precise CKM values are im- 
portant for finding the origin of CP violation 
and revealing new physics beyond the Standard 
Model. There are many avenues to the CKM 
matrix elements open to lattice studies Fig- 
ure |l| shows the Wolfenstein parameterization of 
the CKM matrix with some decays studied on the 
lattice. B semileptonic decays are of particular 
interest as ways of determining \Vub\ and |Vcb|. 



*PIenary talk presented at Lattice 
Australia. 
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Figure 1. The CKM matrix in the Wolfenstein pa- 
rameterization. Decay rate measurements yield 
CKM matrix elements provided hadronic matrix 
elements are known. 



Figure 2. The (p, rf) plane. The semicircular al- 
lowed region centered at (0, 0) (dotted lines) is 
from inclusive B semileptonic decays. The region 
between dashed lines is allowed by B^-B mix- 
ing. The region between solid lines comes from 
CP-violating K^-K mixing. Rosnerj2j. 



The values of \Vub\ and \ Vcb\ from semileptonic 
decays place important constraints upon Wolfen- 
stein parameters p and 77. From Figure ^ the ra- 
tio of CKM elements is A\/ + rf where A is the 
sine of the Cabbibo angle. The ratio scaled by 
A then defines a semi-circle centered at (0, 0) in 
the (/o, ry) plane. Uncertainties in the values of the 
CKM matrix elements widen the semicircle into 
the allowed region shown between the dotted arcs 
in Figure ||. 

The figure also shows allowed regions in the 
(p, 77) plane from B^-B mixing and K^-K mix- 
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ing. With more precise measurements of the 
CKM matrix elements the allowed regions will 
shrink. If all the regions no longer intersect this 
would indicate an inconsistency with the Stan- 
dard Model. Uncertainties in hadronic matrix 
elements are now one of the largest sources of 
uncertainty for CKM matrix elements. 

The ratio |X^6|/|Vcfc| is now best determined 
by looking at the lepton energy spectrum for in- 
clusive decays in the end-point region above the 
charm production limit [^|^. This ratio depends 
upon theoretical form factors. Using current form 
factor models, the Particle Data Group|^ quotes 
IKbl/IKhl = 0.08 ± 0.02 where the 25% error is 
due to combined statistical and theoretical un- 
certainties. Three-family unitarity of the CKM 
matrix and constraints from other CKM element 
measurements yield \Vub\ values that range from 
0.002 to 0.005 at the 90% confidence level. Thus, 
\Vub\ is only known to within a factor of two. 

The CLEO II collaboration has recently re- 
ported the first exclusive measurements for B 
light-meson semileptonic decay modes These 
measurements promise better methods of deter- 
mining \Vub\- Exclusive branching ratios deter- 
minations also require form factors as theoretical 
inputs. With current models the CLEO results 
have an uncomfortably large theoretical uncer- 
tainty. Hence, there is an essential need for re- 
liable hadronic matrix elements. Given the phe- 
nomenological importance of Vub, light-meson de- 
cays modes for the B meson deserve high priority 
in lattice studies. 

The motivation for using the lattice is clear: 
it provides a direct numerical solution of QCD. 
At present, the lattice is the only systematically 
improvable, nonperturbative calculation method 
for QCD. Lattice technology for semileptonic de- 
cays is already well established]^. Systematic er- 
rors must be understood to obtain reliable results 
from the lattice. Lattice discreteness is a source of 
systematic errors. Discretization errors for heavy 
quarks and momentum dependent errors are two 
important potential sources of error considered in 
Section 

Full 0{a) improvement for on-shell matrix 
elements with arbitrary-mass four-component 
quarks is discussed. Improved matrix elements 



for bottom quarks can be studied on present lat- 
tices. Improvement increases the reliability of 
lattice calculations and allows coarser lattices to 
be used. Computations then become faster and 
cheaper. Less costly numerical calculations make 
it cost effective to do systematic studies of other 
sources of error. 

Section discusses momentum dependent er- 
rors. Estimates for momentum dependent errors 
in 0(a)-improved matrix elements on present lat- 
tices show that errors are probably less than 20% 
for momenta below 1 GeV/c. Reducing these er- 
rors would require smaller lattice spacings, adding 
0{a^) improvements to the quark action or ex- 
trapolations to zero lattice spacing at fixed phys- 
ical momentum. 

Section ^ is a survey of D decays on the lat- 
tice. Maximum meson recoil momenta in the D 
rest frame are below 1 GeV/c for semileptonic de- 
cays. Momentum dependent errors are then prob- 
ably below 20% in present studies. Hence, form 
factors at maximum recoil momentum (g^ — 0) 
can be extracted. A summary of lattice D ^ K 
and D —f K* form factors however, shows no 
trend with lattice spacing. Therefore, high preci- 
sion tests for lattice spacing dependence are still 
needed. 

Light-meson decay modes of the B are consid- 
ered in Section ^. On present lattices, one ap- 
proach to B decays is to study D decays and then 
to extrapolate, guided by heavy quark symme- 
try, to the B. Heavy-quark scaling of hadronic 
matrix elements is tested by calculating 0{a)- 
improved matrix elements for charm, bottom and 
static quarks. Matrix elements vary smoothly 
with quark mass. Therefore, when heavy quark 
symmetry is used to find B matrix elements, in- 
terpolations including static matrix elements are 
apt to be more reliable than extrapolations based 
on charm results alone. 

Light mesons have a maximum recoil momen- 
tum of 2.6GeV/c in B decays. Present B cal- 
culations are then restricted to only a fraction of 
the full momentum range. Predictions at large re- 
coil momenta are then potentially unreliable. Dif- 
ferential decay rates and partial widths are sug- 
gested as useful quantities that can be reliably 
calculated with present lattices. 
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Section ^ reports upon decays with heavy 
hadron final states. Studies have focused upon 
mesons in order to understand B D decays and 
extract \Vcb\- Baryon decays are now also under 
study. Baryons offer a new laboratory for heavy 
quark symmetry studies - including baryon Isgur- 
Wise universal form factors. Finally, two studies 
of Heavy Quark Effective Theory on the lattice 
are discussed. One study examines nonpertur- 
bative renormalization for lattice Heavy Quark 
Effective Theory while the second demonstrates 
a dramatic reduction in statistical noise when us- 
ing optimized wavefunctions in calculations of the 
meson Isgur-Wise function. 

2. ARBITRARY-MASS QUARKS 

A quark is considered heavy when mg ^ ^qcd 
as is the case for charm and bottom quarks. 
Wilson and Sheikholeslami-Wohlert (SW) heavy 
quarks are widely used in matrix element stud- 
ies. On present lattices heavy-quark masses are 
typically a large fraction of the lattice cutoff 
Too ^ 1/ a. Errors arising from the discrete lattice 
can then be important. Standard Wilson quarks 
have 0(a) errors in matrix elements. Since de- 
creasing the lattice spacing so that utoq <C 1 for 
a bottom quark is extremely costly, improvements 
which remove discretization errors on present lat- 
tices become an attractive alternative. This sec- 
tion focuses on a method for systematically re- 
moving lattice-spacing errors in four-component 
quarks of arbitrary mass. 

Lattice-spacing errors can be observed in re- 
sults for pseudoscalar decay constants. Figure |^ 
shows decay constants from a study using both 
Wilson and SW quarks on a lattice where the D- 
meson mass is amjj ^ Wilson results are 

indicated by the diamond symbols in the figure. 
Conventional Wilson quarks incorrectly lead to 
af^/am = in the infinite-mass limit. Hence 
discretization errors overwhelm Wilson quarks at 
large masses. The figure shows that extrapola- 
tions of Wilson results from the D to the B are 
apt to be unreliable. 

Errors in on-shell matrix elements can be re- 
duced systematically by improving the lattice ac- 
tion and currents. 0{a) (Clover) improvement 
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The quantity afp^anip for heavy 
Three types of quarks are com- 
pared: conventional Wilson quarks (diamonds), 
improved Wilson quarks (octagons) and Clover 
quarks (squares). At large mass conventional 
Wilson quarks have large discretization errors. 
The D meson mass is {amD)~^ « 1 on this lat- 
tice. 



for light quarks is based upon the SW action[gj. 
The square symbols in Figure ^ denote Clover 
decay constants. Clover results extrapolate to 
an infinite-mass limit closer to the static value 
afy^am « 0.32 indicating a sizable reduction in 
discretization errors at large mass compared to 
Wilson quarks. 

Consider again the Wilson results. Large 
discretization errors are apparently removed by 
adopting an unconventional wavefunction renor- 
malization for Wilson quarks 1 10, Results for 
reinterpreted Wilson quarks (circles in Figure H) 
also have an infinite-mass limit in better agree- 
ment with the static result. With the uncon- 
ventional wavefunction renormalization the cor- 
rect free-field quark propagator is recovered at 
infinite-mass. Since improvement is possible in 
both heavy and light limits, the Fermilab group 
then investigated the possibility of improvements 
for quarks of arbitrary mass. 

The Fermilab program improves arbitrary- 
mass four-component quarks |l2]. It is equiva- 



4 



lent to the usual improvements for light quarks 
and to the static approximation at infinite mass. 
Fully 0(a)-improved quarks are obtained using 
the SW action and improved currents. Bottom 
quarks can be put on existing lattices without 
0{a) errors. Heavy Wilson and SW quarks are 
used successfully to study -onia[^. In B-decay 
studies uncertainties due to extrapolations in the 
heavy-quark mass can be eliminated and calcu- 
lations become cheaper and faster since present 
lattices can be used. Computational savings can 
then be used to remove remaining discretization 
errors by repeating calculations at several lattice 
spacings. 

2.1. The Action 

The 0{a) quark action in the Fermilab formal- 
ism is a sum of three terms: 

Sf^So + SB + Se. (1) 

Term 5*0 has dimension-three and -four operators 
and a dimension-five Wilson operator to eliminate 
doublers: 



(2) 



Operators Z?j^ are the forward (+) and backward 
(— ) first differences, is the central difference 
and As is the spatial laplacian. By setting — 
1 and C = 1 action Sq becomes the Wilson action. 
The action is asymmetric in space and time when 

The chromomagnetic and chromoelectric terms 
are: 

Sb = -^acsC y"*'S-B*' (3) 

Se = -\acEC (4) 

In contrast to the light-quark action where cb — 
Ce and Sb + Se is "clover" term, here cb and ce 
may be unequal. 

Parameters Vs Ci cb and ce are functions of 
the coupling; tree-level values are known. Higher- 
order corrections can be obtained perturbatively. 



These parameters are also functions of quark 
mass. The mass dependence is treated to all or- 
ders. 

Values = 1, C = 1 and c ~ cb = ce are used 
for light quarks. With c — 1 the action is the tree- 
level SW action; tadpole improvement increases 
c to between 1.4 and 1.5 on typical lattices. The 
tadpole-improved SW action leads to fully 0{a)- 
improved light quarks. 0{a) improvement is still 
obtained for heavy quarks provided currents and 
the physical mass are suitably defined. 

2.2. Physical Heavy-Quark Masses 

In the non-relativistic limit, the Hamiltonian 
corresponding to Sf is 



H 



* Ml + 7oylo - 



2M2 



• B 

2Mb 



(5) 



where 4* is the physical field defined in Equa- 
tion H. This looks like the Pauli Hamiltonian; 
Ml, M2, Mb, and Me however are functions of 
toq. This Hamiltonian describes a non-relativistic 
quark when M2 = Mb = Me- For both Wil- 
son and SW quarks, when amg 1 the kinetic 
mass, M2, and the rest mass. Mi, are signifi- 
cantly different. Since Mi is unimportant for non- 
relativistic quarks, M2 should then be considered 
the physical mass. 

The Hamiltonian is a useful guide for estimat- 
ing errors for non-relativistic quarks. Systematic 
errors will arise if Mb and Me do not equal M2. 
At tree-level, when amo ^ 1 Wilson quarks have 
Mb > M2 while SW quarks have Mb ^ M2. 
Therefore, from Equation ^, Wilson quarks can 
be expected to correctly reproduce only spin- 
averaged energy levels while SW quarks should 
also approximate spin-dependent features. 

Physical quark masses can be found non- 
perturbatively in -onia by using the non- 
relativistic energy dispersion relation 



-2 

^1 -mi + - — 



(6) 



The meson kinetic mass, TO2, and the rest mass, 
TOi, differ significantly when aniQ 1 for Wil- 
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Figure 4. Meson kinetic mass, m2, as a function 
of the rest mass, mi in -onia. Points are lattice 
determinations using the energy dispersion rela- 
tion. The curve is based upon the tree-level ex- 
pression for the quark kinetic mass in terms of 
the rest mass. 



son and SW quarks. Charm and bottom quark 
masses are determined by matching the physical 
mass and m2- Spin-averaged levels are used; the 
physical lattice spacing may be found using the 
IP- IS* mass splitting. 

Figure ^ shows am2 and ami values obtained 
for three lattice spacings using SW quarks. Nu- 
merical results agree well with the curve Q which 
is based on the tree-level functional relation be- 
tween quark masses M2 and Mi. Differences 
between numerical results and tree-level predic- 
tions are, in part, due to higher-order corrections. 
One-loop corrections to quark masses are being 
computed Ip^. The figure shows that the rest 
mass systematically underestimates the physical 
mass. At /3 = 6.1, M2/M1 is about 1.2 for the 
charm quark and 2.8 for the bottom quark. Us- 
ing the rest mass to determine the physical mass 
rather than the more physical kinetic mass intro- 
duces a systematic error into heavy-quark calcu- 
lations. 

2.3. Currents 

Fully 0(a)-improved matrix elements require 
0(a)-improved currents. Physical fermion opera- 
tors, ^, appear in renormalized currents. These 



Tree — Level dj (Tadpole— Improved) 
= 5.9 




Figure 5. Tadpole- improved tree-level values of 
coefficient di plotted as a function of I/A/2 for 
/3 — 5.9. Coefficient di equals zero in the infinite- 
and zero-mass limits and is typically largest be- 
tween charm and bottom. 



fields are constructed from fields ^P", which ap- 
pear in the hopping parameter form of action Sf. 
Fields 4" are given by lira 



1 -I- adi(mo, 5^) 7 • D -t- 0(a2) 



* = V2Kexp(aA/i/2)^'. 



(7) 
(8) 



The three-gradient term is an 0{a) correction to 
the current; the exponential factor is a correction 
in all-orders of arriQ . 

Coefficient di depends upon the action and is 
a function of the coupling and quark mass. Fig- 
ure [s] shows di as a function of 1 /M2 for tadpole- 
improved SW quarks at /3 = 5.9. At tree-level, 
di vanishes in both zero- and infinite-mass lim- 
its. The correction is typically largest between 
charm and bottom. At f3 = 5.9, lattice results 
show that /d increases by about 7% when the 
three- gradient term is included in the current. 

In Equation H the factor exp(aMi/2) appears in 
addition to the conventional \/2k factor for Wil- 
son and SW quarks. This factor removes large 
discretization errors when aniQ ^1. Its effect 
upon Wilson decay constants is seen in Figure ||. 
To O(amo) this factor is equivalent, at tree- level, 
to the Clover rotation used to improve light- 
quark currents. The factor becomes vl — 6k us- 
ing the tree-level expression for the rest mass and 
•\/l — 6uqk with tadpole improvement. The aver- 
age plaquette or kappa critical is typically used 
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Figure 6. All-order corrections beyond tree-level 
for the vector current renormalization. Correc- 
tions have a mild dependence upon quark mass. 
They are consistent with perturbation theory for 
small masses. At infinite mass higher-order cor- 
rections are zero and the tadpole-improved tree- 
level renormalization yields the exact result. 



to estimate uo|17|]. 

An 0(a)-improved current coupling quarks Q 
and q has the form 



(9) 



where the fields are given in Equations ^ and |§. 
Matrix F'^ equals 7^ for a vector current and 7'^ 7^ 
for an axial current. Renormalization Zj matches 
the lattice current to the continuum. 
With tadpole improvement, Zj becomes 



Zj = ^^kqlfiq exp 



aMl) 



(10) 



where Z^ is the vertex correction. Here, fc = mqk 
and Ml denotes tadpole-improved M\. Zj keeps 
this form to any order in . For example, at 



0{g ) expansions Mi = M 

.2 



g^Mf ' and Zr = 
The higher- 



1 + g Zf, appear in Equation [10, 
order corrections to Zj are also functions of the 
quark masses. 

The mass dependence of higher-order terms in 
Zj can be checked non-perturbatively for the vec- 
tor current by examining the charge normaliza- 
tion of three-point functions. Figure |^ shows 
the ratio of the tree-level renormalization fac- 
tor, Zy'^'^, and the nonperturbative factor, Zy'', 



versus I/M2 for tadpole-improved SW quarks at 
(3 ~ 5.9. Terms in Zy^ beyond tree level cause 
this ratio to differ from unity. The figure shows 
that higher-order corrections change Zy^'^ by at 
most 20% and that light-quark corrections are 
largest. One-loop corrections of around 20% are 
expected for light quarks at /3 = 5.9. Note 
that renormalization Zy^^ is exact in the static 
limit [T^ . Thus, higher-order corrections depend 
smoothly on mass and yield the expected results 
in the light- and infinite-mass limits. 

The Wuppertal Group ||l^ has also checked the 
renormalization in 
Equation |lO| non-perturbatively. They look at ra- 
tios of the local-current three-point function over 
the conserved-current three-point function. They 
test Equation |l^ as function of quark mass using 
flavor-conserving currents, heavy-light currents, 
and currents carrying momentum. The last two 
cases are important tests of the current renormal- 
ization for semileptonic decays. 

They find that Zy'^'^, given by Equation |l^, 
is in reasonable agreement with the nonpertur- 
bative determination of Zy for the three cases 
studied. Because of mass-dependent higher-order 
terms, exact agreement between the nonperturba- 
tive renormalization and Zy*^*^ is not expected. In 



contrast to Equation 10, mass- independent renor- 
malization using the perturbative result for light 
quarks leads to poor agreement with the nonper- 
turbative results for heavy quarks. 

2.4. Momentum Dependent Errors 

With fully 0(a)-improved actions and currents 
0{ap) errors are removed from matrix elements. 
Since improvement is carried out for on-shell 
quarks at rest, however, momentum dependent 
discretization errors are still a concern at large 
momenta \ap\ 3> 0. It is then important to con- 
sider the effect of momentum dependent errors 
in semileptonic decays. Tree-level errors in quark 
matrix elements can be determined analytically 
by comparing infinite-volume lattice matrix ele- 
ments and continuum matrix elements. Errors 
in the quark matrix elements can then serve as 
a guide to estimating errors in hadronic matrix 
elements. 

Consider matrix elements for an external vector 
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Figure 7. Momentum dependent errors in tree-level quark matrix elements (c(p3;, t)|v^|6(o,T)) for typical 
lattice spacings. Errors are shown as a function of the charm-quark physical momentum. Plots show 
errors in temporal (top) and spatial (bottom) matrix elements. Currents are the local current (left) and 
the 0(a)-improved current (right). 



current {q{pf,crf)\V^\Q{pi,ai)) where the cur- 
rent is defined in Equation y. Matrix elements for 
the local current {di = in Equation 0) and for 
the 0(a)-improved current are calculated. Wilson 
quarks with equal kinetic and physical masses are 
used. For simplicity, the initial quark is at rest. 
This should be a good approximation to a heavy 
quark decaying within a meson at rest. Since ma- 
trix elements are independent of the mass of the 
at-rest quark, error estimates are the same for 
charm and bottom decays. Here, the initial quark 
is taken to be a bottom quark. Only on-axis recoil 
momenta are considered for the final quark. 

shows the relative error in lattice 



Figure 2.4 



matrix elements with a charm-quark final state 



as a function of the physical momentum. Plots 
on the left show errors for the local current, 
while those on the right show errors for the 0(a)- 
improved current. Upper plots correspond to 
temporal matrix elements while the lower plots 
show spatial matrix elements. Plots show typical 
lattice spacings in the range 0.17 to 0.06 fm. 

Errors change more systematically with mo- 
mentum with the improved current. Hence dis- 
cretization errors in hadronic matrix elements are 
likely to be under better control at small mo- 
menta with the improved current. Note however 
that even with full 0(a) improvement errors in 
spatial matrix elements may still be of order 20% 
at \p\ = 1.24GeV/c on a /? = 5.9 lattice. A mo- 
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mentum of 1.24GeV/c corresponds to two units 
of lattice momentum on a two fermi lattice. Er- 
rors can be reduced by adding O(a^) improve- 
ments to the action and currents. Note however 
that precision matrix elements may still be ob- 
tainable in the presence of sizable discretization 
errors by extrapolating lattice matrix elements to 
zero lattice spacing at fixed physical momenta. 
The feasibility of this approach hinges upon us- 
ing suitable techniques to control statistical noise 
in lattice correlators with momentum pO|] . 

3. D DECAYS 

Semileptonic decays on the lattice were pio- 
neered with studies of D mesons ]2l| , p2| . D decays 
continue to be an active area of investigation! 23 
|2|,|25| Mp8| . At this conference the Los 
Alamos[2^and Wuppertal[p^ groups presented 
updates for their D studies. Charm decays are 
an interesting laboratory for QCD and serve as 
a first step towards the study of B decays. D- 
decay experiments provide a basis for comparison 
in lattice studies. Experimental checks for D de- 
cays are an important way of checking systematic 
errors in lattice calculations. 

Experimental information for D decays has 
grown steadily since the first lattice calculations. 
Nevertheless, the K and K* recoil momentum 
distributions are still not well known p9[. Instead, 
pole dominance is typically used to extract form 
factor values at maximum recoil momentum cor- 
responding to = [pD — Pk)^ = 0. For the K, 
the location of the pole can also be determined 
from the data. Experiments find pole masses 
consistent with identifying the D* as the dom- 
inant pole. For K*, poles masses are fixed to 
the D* and D** masses and form factor ratios 
Ry = V{0)/Ai{0) and i?2 = A2(0)/Ai(0) are 
found. Individual vector form factors are then 
found by assuming the branching ratio depends 
predominantly on A^. Hence experimental evi- 
dence for pole dominance in D decays is indirect. 

Lattice calculations for D decays cover the 
whole range of recoil momentum. Statistical and 
systematic errors at maximum recoil are impor- 
tant to estimate since D results are tradition- 
ally presented as form factors at maximum re- 



coil. In the D rest frame the maximum kaon mo- 
mentum is 0.86GeV/c while the maximum pion 
momentum is 0.93GeV/c. Statistical errors at 
such momenta can be minimized by using good 
source operators for moving mesons. With ~ 300 
gauge configurations statistical errors below 10% 
are likely for the whole recoil-momentum range. 

Momentum dependent errors at maximum re- 
coil can be estimated from errors in tree-level 
quark matrix elements. Errors in fi^{0) are pre- 
dicted to be around 10% with 0(a)-improved 
quarks at (3 = 5.9. This and other discretization 
errors can be reduced by working on a finer lat- 
tice. With improved actions, however, extremely 
fine lattice spacings are not necessary. A series 
of faster less costly computations on relatively 
coarse lattice spacings can be done. Remaining 
discretization errors are then removed by extrap- 
olation to zero lattice spacing pO||. 

Table |l] is a summary of D form factors at 

= obtained in lattice studies. Statistical er- 
rors are shown with estimates, when given, for 
systematic uncertainties. Best estimates of form 
factors from experimental studies are shown for 
comparison To investigate lattice-spacing er- 
rors, lattice spacings, lengths and the quark ac- 
tion used in each study are listed in the table. 
Lattice results are mostly within two sigma of the 
experimental values. Unfortunately, while statis- 
tical errors for are typically below 15% sta- 
tistical errors for the other form factors tend to 
be larger. Hence, from the tabulated results it is 
difhcult to determine how lattice-spacing errors 
affect form factors. Higher precision studies of D 
decays are needed. Comparison with experiment 
will be an important test of whether systematic 
errors are under control. 

Since discrete lattice momenta typically 
bracket q^ = 0, the form factors in Table |l] are ob- 
tained by interpolation or, in cases where statis- 
tical errors were too large, by extrapolation over 
a limited range of recoil momentum. In prac- 
tice, pole dominance provides a simple means of 
smoothly interpolating lattice results. 

The Los Alamos and Wuppertaljl^ groups 
both test possible interpolation functions - in- 
cluding pole dominance forms. The Los Alamos 
group note that at leading order in l/rric the 
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Table 1 

D ^ Klv and D K* Iv form factors at = 0. Experimental averages (EX) and a summary of lattice 
results (LA) are shown. Lattice spacings, a, lengths, L (both in fm) and the quark action (Ac) used in 
each study are also shown. 



Reference 



Ac 



f+(0) 



V(0) 



Ai(0) 



EX 



Averagej 29| 
LMM S M 
BKS 124 1 
ELC 1251 
UKQCD 
LANL H? 
LAT-APE 
W'taliTl 



0.75 ±0.03 



1.1±0.2 



0.56 ±0.04 



LA 




.095 
.095 
.055 
.071 
.095 
.095 
.064 



1.0 
2.3 
1.3 
1.7 
3.0 
1.7 
1.5 



W 
W 

w 

sw 

w 

sw 

w 



0.63 ± 
0.90 ± 
0.60 ± 
0.67i;] 
0.71 ± 
0.78 ± 
0.71 ± 



0.08 

0.08 ±0.21 
0.15 ±0.07 

,07 
.08 

0.04 
0.08 



0.12 



+0.10 
-0.07 



0.86 ±0.10 

1.43 ±0.45 ±0.49 

0.86 ±0.24 

-■^■UJ^-O.lS 

1.28 ±0.07 
1.08 ±0.22 
1.34 ±0.24 



+0.19 
-0.14 



0.53 ±0.03 

0.83 ±0.14 ±0.28 

0.64 ±0.16 

'^-0.10 

0.72 ±0.03 
0.67±0.11 
0.61 ±0.06 



0.09 
-0.07 



Reference 



Ac 



A2(0) 



V(0)/Ai(0) 



A2(0)/Ai(0) 



Average ^ 



EX 



0.40 ±0.08 



1.89 ±0.25 



0.73 ±0.15 



LA 



LMMS |2§ 
BKS |24| 
ELC ll 
UKQC D 
LANL 1^ 
LAT-APE H 
W'taliTl 



095 


1.0 


W 


095 


2.3 


W 


055 


1.3 


w 


071 


1.7 


sw 


095 


3.0 


w 


095 


1.7 


sw 


064 


1.5 


w 



0.19 ± 
0.59 ± 
0.40 ± 
0.66l^ 
0.49 ± 
0.49 ± 
0.83 ± 



0.21 

0.14 ±0.24 
0.28 ±0.04 

10 
.15 

0.09 
0.34 



1.6±0.2 

1.99 ±0.22 ±0.33 
1.3 ±0.2 

1 4+0-5 

1.78 ±0.07 
1.6 ±0.3 



0.4 ±0.4 

0.7 ±0.16 ±0.17 

0.6 ±0.3 

0.9 ±0.2 

0.68 ±0.11 

0.7 ±0.4 



0.20 



+0.12 
-0.08 



heavy quark symmetry relation between D form 
factors /+ and /o appears to be incompatible with 
both form factors having a pole form. They find 
favors a pole mass lighter than the D* mass 
while /o yields the expected pole mass. To de- 
termine /:f^(0) they use the pole form giving the 
lighter effective pole rather than fixing the pole 
mass to the D* mass. 

Although functions inspired by pole dominance 
are commonly used for interpolations of lattice D 
form factors it is important to remember that how 
adequately pole dominance describes D decays or 
especially B decays is not known. 



4. B DECAYS TO LIGHT MESONS 

The light-meson decay modes for the B are es- 
pecially important calculations for the lattice now 
that the CLEG II collaboration has branching ra- 
tios for the pion and rho|^. With a sample of 2.6 
million B B decays they find branching ratios 



(scaled by x 10 ^) 



6(5° 
S(B0 



l+v) 



1.34 ±0.35 ±0.28 (ISGW) 
1.63 ±0.46 ±0.34 (BSW) 

p-l+v) - 

2.28 ± 0.36 ± 0.59 (t;^;!^^) (ISGW) 
3.88 ± 0.54 ± 1.01 (t[3:?^) (BSW)' 

The first error for each ratio is statistical and 
the second is systematic. The third error for 
the rho is the uncertainty due to non-resonant 
TTTT contributions. These branching ratios depend 
upon theoretical input for the efficiency calcu- 
lation. The branching ratios above result when 
the lSGW|l| and the WSB|3|] quark models 
are used for efficiencies. Differences between the 
models indicate that theoretical uncertainties in 
the hadronic matrix elements are a substantial 
source of error. Note that reliable hadronic ma- 
trix elements are even more important for deter- 
mining |Kb|- 

Light decay modes are under investigation on 
the lattice p^ , ^ . ^ 34 1^. These studies all use 
either standard Wilson or Clover]^ quarks. Be- 
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cause of concerns over discretization errors, direct 
calculations involve heavy mesons with masses 
around the D mass. Heavy quark symmetry is 
then used to guide extrapolations from the D sys- 
tem to the B meson. 

Hadronic matrix elements for light-meson de- 
cay modes scale in the heavy-quark limit [ ^J36| . 
For example, consider decays H ^ nlv. In the 
i?-meson rest frame 



[rriH 



,2//3o 



it{p^)\V^^\H{Q) 



M{p.) (11) 



where M. is independent of the heavy-meson 
mass. The relation is not an equality since there 
are 0{l/mH) corrections. In particular, the recoil 
momentum is assumed to be small \p-7:\/rnH ^ 1. 
Equation |ll| relates D and B matrix elements 
with equal pion momentum. If 0{l/mH) scale- 
breaking terms were small this expression would 
provide useful information about B decays from 
experimental data for D decays. 

Form factor scaling relations follow from ma- 
trix element relations. For example, from Equa- 
tion |l| one finds /+ / y^mn ^ + S+/ mn ) 
(modulo logarithms) where the linear l/mn 
scale-breaking term is now explicit. On the lat- 
tice, linear functions are typically used for extrap- 
olations. In a pioneering study using this method 
the ELC group found 1/mjj terms which differ 
from zero by at most two sigma[^. Statisti- 
cal errors however did not exclude slopes as large 
as 1 GeV for /+ / y/niH and 0.7 GeV for Ai^ttih- 
These statistical errors and possible systematic 
errors in heavy-quark extrapolations are sources 
of uncertainty in S-meson results. 

The quantity (tt jF'^l H) /^/2Mh 2E^ in Equa- 
tion ^ is obtained directly in lattice calculations 
when the heavy meson is at rest. The depen- 
dence upon heavy-quark mass has been investi- 
gated by the Fermilab group on a /3 = 5.9 lattice. 
Matrix elements were computed using fully 0{a)- 
improved charm and bottom quarks with strange- 
mass light quarks. Results were also obtained 
in the static approximation. Figure || shows the 
hadronic matrix elements parameterized by pion 
momentum and plotted versus l/mn- Plot sym- 
bols and dashed curves distinguish pion momen- 
tum values. Momenta in lattice units are indi- 



cated in the plots; one unit of momentum cor- 
responds to 0.70 GeV/c. The plot on the left 
shows temporal matrix elements and on the right 
are x-component matrix elements. Matrix ele- 
ments vary smoothly from charm to the static 
limit. Therefore in studies employing extrapola- 
tions from the D to the B, results are likely to 
be more reliable when static matrix elements are 
included so that B results are then obtained by 
interpolation. 

Light-meson decay modes for the B have a 
maximum recoil momentum of 2.6 GeV/c com- 
pared to only 0.93 GeV/c for D decays. Concerns 
over momentum dependent discretization errors 
and statistical errors which increase rapidly with 
momentum limit current calculations to momenta 
around 1 GeV/c. Hence in contrast to D decays, 
B form factors at 5^ = are far beyond the range 
of direct calculation. In addition, extrapolations 
to large recoil momenta are also likely to be unre- 
liable. Extrapolations to large momenta can am- 
plify small discretization errors which occur for 
small momenta. Also, extrapolations which as- 
sume pole dominance for the B introduce an un- 
known degree of model dependence. Therefore it 
is best to focus upon quantities in B decays that 
can be directly calculated without extrapolations. 

Differential decay rates and partial widths do 
not require large extrapolations in momentum. 
Consider the differential decay rate for B ^ irlv 



dV 



= 2ms 



247r3 



(12) 



d\p^\ 24^3 

Figure ^ shows this differential rate (less the mo- 
mentum independent pre-factors) in a prelimi- 
nary study of B decays by the Fermilab group. 
Statistical errors are shown by the solid-line error 
bars; larger (broken-line) error bars show momen- 
tum dependent errors estimated using tree-level 
quark matrix elements. In this /3 = 5.9 study 
statistical errors are below 10% and momentum 
dependent errors are estimated to be less than 
15% below 1 GeV/c. Errors may then be small 
enough that extrapolations to zero lattice spacing 
are reliable enough to remove most discretization 
errors. 

Results in Figure ^ are for strange-mass light 
quarks. Chiral extrapolations are necessary to 



11 



1.5 



1.0 



0.5 
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\ \ \ \ \ \ 



0.0 



4 (0, 0, 0) 



-4- 



t 



0.5 
1 / 
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1 / 
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(1. 0) 
T (2, 0. 0)_ 

--4 (1, 1, 1) 



i 

m„ 
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Figure 8. Hadronic matrix elements for static, bottom, and charm quarks at /3 = 5.9. The plot on the left 
shows temporal matrix elements. On the right are x-direction matrix elements which are proportional 
to the x-component of pV- Curves and distinct plot symbols indicate pion momenta which are shown in 
lattice units. 



1.0 



> 



0.5 



dr(B, -> lu) 



d|pj 



0.0 





1 



_ 2 
IpJ (GeV/c) 
Figure 9. B-meson differential decay width for 
the light-pseudoscalar mode. Light-quark masses 
are approximately equal to the strange quark 
mass. Solid error bars are statistical errors. Es- 
timated momentum dependent errors are shown 
by the larger broken- line error bars. 
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obtain decay rates for the physical pion. Chi- 
ral extrapolations of form factor /+ at small pion 
momenta may be difficult since soft pion argu- 
ments suggest /+ varies rapidly with pion mass in 
heavy- meson decays js^. Note however variations 
in /+ at small recoil momenta are suppressed by 
four powers of the momentum in the decay rate. 
Hence, chiral extrapolations of differential decay 
rates may be more robust than chiral extrapola- 
tions of form factors. 

Partial widths can be obtained by smoothly 
interpolating between lattice differential decay 
rates values and then integrating. Partial widths 
obtained this way are expected to be relatively 
model independent for momenta within the range 
of direct lattice calculations. The UKQCD col- 
laboration has found widths using this method 
for the vector decay mode[Q. Table || shows 
their results for differential decay rates and par- 
tial widths along with statistical errors. Chiral 
extrapolations have not been done and strange- 
mass light quarks are used. Note that at 
Qq = 14.4(GeV/c)^ which corresponds to \pp\ = 
1.2 GeV/c statistical errors are about 10% for the 
partial width. The authors point out that partial 
widths for the rho and pion decay modes from 
the lattice can be used to extract \Vub\ in a model 



12 



Table 2 

B-meson differential rates and partial widths for 
the light-vector decay mode. Momenta are in 
GeV/c. Width T is in units of iKbP x IQ-^^ GeV. 
Strange-mass light quarks are used. UKQCD 
collaborationltMl . 





dT/dq%2 




20.3 


0.00 


0.0 


19.7 (t}) 


0.19(1^) 


0.08 {+{) 


17.5 (tl) 
16.7 (tl) 

15.3 HI) 

14.4 (11) 


0.57 (+1) 
0.58 {tl) 
0.6(11) 
0.8(1?) 


0.9(11) 
1.3(11) 
2.3 (1^) 
3.0 (1^) 


charm threshold 




11.6 
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independent way from exclusive B decay modes. 

The CLEO II exclusive- mode analysis suggests 
that a sizable fraction of pions have momenta be- 
low 1.4GeV/c. Hence exclusive partial widths 
from the lattice may in fact be a viable means of 
extracting \Vub\- Minimizing errors in \Vub\ will 
require adjusting the maximum recoil momentum 
limit in the partial width so as to minimize both 
lattice and experimental errors. 

The most important task for semileptonic de- 
cay studies is improving the precision to which 
\Vcb\ is known. It is then important to eliminate 
potential sources of error. Extrapolations to large 
recoil momenta using present lattice results are 
unreliable. Differential decay rates and partial 
widths are likely to be more reliable and present 
lattice results in a form more useful to the exper- 
imental community. 

5. DECAYS TO HEAVY HADRONS 

Heavy-meson decay modes are studied to test 
non-perturbatively heavy quark symmetry and to 
provide theoretical input needed to extract \Vcb\ 
from _B-meson decay rates. The study of these 
decay modes on the lattice is discussed in detail 
by Kenway[^ and by Lellouch|Q. This section 
provides a update in the form of results presented 



at this conference. 

The Isgur-Wise function (up to 0(l/m) cor- 
rections) has been extracted from QCD ma- 
trix elements calculated using Wilson and SW 
quarks 1^,^. Systematic error estimates for 
these calculations are not yet complete. Tree- 
level quark matrix elements however provide a 
way of estimating momentum-dependent errors. 
Note that according to Figure 2.4 relative errors 
are smaller for temporal matrix elements. Hence 
elastic form factors which require only tempo- 
ral matrix elements are likely to yield the small- 
est momentum-dependent uncertainties for the 
Isgur-Wise function. The figure shows that tree- 
level errors are probably less than 10% for mo- 
menta below IGeV/c. It would be interesting 
to see how errors bias the determination of the 
Isgur-Wise slope at the zero-recoil point. 

Another means of studying heavy quark sym- 
metry is to formulate Heavy Quark Effective 
Theory (HQET) directly on the lattice @. At 
leading order lattice HQET (LHQET) describes 
an infinitely massive quark with arbitrary four- 
velocity. The lattice theory must be renormalized 
to match HQET. 

At this conference, Mandula and Ogilvie in- 
vestigate velocity renormalizationji^. In HQET 
the multiplicative renormalization of the heavy 
quark's four-velocity is actually a wavefunc- 
tion renormalization. Lorentz symmetry is re- 
duced on the lattice to hypercubic symmetry. 
Then, in LHQET additive shifts can also oc- 
cur. These renormalizations have been computed 
perturbativelylQ. The perturbative expecta- 
tions are compared with nonperturbative deter- 
minations of the renormalization coefficients from 
the lattice. In this pilot study, the coefficient of 
the term proportional to four-velocity is within 
one sigma of the one-loop value. Perturbative 
and lattice coefficients for terms cubic in the four- 
velocity were in complete disagreement however. 
Tadpole improvement of the LHQET action and 
a Lepage-Mackenzie treatment of perturbation 
theory |l 711 may help resolve these discrepancies. 

The Kentucky group presented preliminary nu- 
merical results for the unrenormalized Isgur-Wise 
form factor in LHQET [Q. They find their varia- 
tional method (MOST) for constructing optimal 
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wavefunctions can be effectively used in LHQET. 
MOST operators liave good signal-to-noise ratios. 
With 32 gauge configurations statistical errors are 
around 5% while errors due to excited state con- 
taminations are estimated to be 5%. These pre- 
liminary results encourage further study. Com- 
parisons of the Isgur-Wise function directly com- 
puted using LHQET and the Isgur-Wise function 
extracted using finite- mass quarks will provide in- 
formation about 0(l/m) errors. 

Baryons containing one heavy quark are under 
study by the UKQCD collaboration!^. Semilep- 
tonic decays of baryons may provide additional 
checks on CKM matrix elements. Baryons also 
offer additional tests of heavy quark symmetry. 
Consider the six form factors for Ab Ac liy. 
They are related to a single baryon Isgur-Wise 
function appearing in the heavy quark limit. The 
UKQCD collaboration presented preliminary re- 
sults for the axial form factor which is expected 
to have small 0(l/m) corrections. They find lit- 
tle dependence upon 1/m when the heavy-quark 
mass is varied about the charm-quark mass. 

Heavy hadron decays are important for deter- 
mining \Vcb\- Lattice calculations provide valu- 
able nonperturbative tests of heavy quark sym- 
metry as well as Isgur-Wise form factors. HQET 
on the lattice is still being explored. In LHQET 
calculations statistical noise and poor ground- 
state isolation are important problems to be over- 
come. Highly optimal particle operators are then 
necessary. 

6. REMARKS 

Semileptonic decays provide important infor- 
mation about CKM matrix elements involving 
heavy quarks. Recent CLEO II results for the 
light-meson decay modes of the B are exciting 
since they promise better determinations of | Vub \ ■ 
Since reliable matrix elements are a crucial to 
this goal, B decays deserve high priority in lattice 
studies. 

Systematic errors must be understood to ob- 
tain reliable results from the lattice. Discretiza- 
tion errors are an important source of uncertainty. 
Improvement docs not require expensive lattices 
with drastically smaller lattice spacings to de- 



crease discretization errors. Improvement is pos- 
sible with arbitrary-mass four-component quarks. 
Bottom quarks can then be put on existing lat- 
tices. Extrapolations in the heavy-quark mass 
which are a potential source of systematic error 
in bottom-quark studies can then be eliminated. 

With 0(a) improvement momentum dependent 
discretization errors in matrix elements are still a 
concern on present lattices for momenta above 

I GeV/c. Alternatives to drastically reducing the 
lattice spacing as a way of controlling momen- 
tum dependent errors include: O(a^) improve- 
ment, and extrapolations to zero lattice spacing 
at fixed physical momenta. 

Present D decay calculations cover the full 
range of recoil momentum. Careful attention to 
sources of systematic error should lead to reliable 
matrix elements. Experimental data for D de- 
cays provide a basis for comparing lattice results. 
Checks for the D are important for testing that 
systematic errors are under control. 

For light-meson decay modes of the B the max- 
imum recoil momentum is greater than the largest 
momentum than can be confidently studied with 
present lattices. Attention should then be di- 
rected towards quantities that can be reliably de- 
termined. Differential decay rates and partial 
widths are a useful form for lattice results that 
do not require extrapolations to large recoil mo- 
menta. Partial widths from the lattice may offer 
a viable means of extracting \ Vub\ from the CLEO 

II exclusive decay measurements. 
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